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Abstract 

Appreciable zero-energy surface bound states are shown to exist in Fe-pnictide superconductors 
with antiphase s^-wave pairing. In contrast to the [110] bound states in d-wave superconductors, 
these bound states arise from the non-conservation of momentum perpendicular to the interface for 
tunneling electrons and the pairing, and hence they can only exist in a small window (~ ib6°)in 
the orientation of edges near [100] direction. Our results explain why zero-bias conductance peak 
was often observed in tunneling spectrum and when it disappears, two coherent peaks show up, 
providing unambiguous signals to test the possible s^-wave pairing in Fe-pnictide superconductors. 

PACS numbers: 74.20.Rp,74.45.+c,74.50.-Hr 
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Recently superconductivity has been observed in several classes of Fe-pnictide materi- 
als [liJsl. One key issue towards understanding their superconductivity lies in the pairing 
symmetry of the Cooper pairs. A conclusive observation of the pairing symmetry remains 
unsettled for iron pnictides to which their superconducting (SC) gap symmetries are ob- 
served to be remarkably dependent of material classes and doping levels j4-6|. Among 
various candidates, the most natural and promising pairing state is considered to be the 
s^-wave which has a sign reversal between a and /3 bands and can be naturally explained 
by the spin fluctuation mechanism jr-ll]. 

Among phase-sensitive measurements, the point-contact Andreev-reflection spectroscopy 
(PCARS) is one of the high-resolution probes for detesting the SC pairing state. The zero- 
bias conductance peak (ZBCP) associated with Andreev bound states (ABS) in PCARS 
measurement has given a direct evidence on the d-wave pairing of high- Tc cuprate super- 



conductors 
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161]. However, PCARS measurements on Fe-based superconductors has yet 



not yielded consistent results. Some PCARS measurements show two coherent peaks and 
indicate that SC pairing state might be fully gaped on the Fermi surface (FS) Other 
measurements show a ZBCP and imply the presence of zero-energy bound states or ABS 
on the interface 
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22| . Moreover, depending on the direction of the sample intersurface. 
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some PCARS has shown ZBCP coexisting with finite-energy coherent peaks 

To address these puzzling observations, surface bound states (or midgap states) have been 
studied and found to appear at both zero and nonzero energies when considering the quasi- 
particles (QPs) interplaying between the two bands 25|-|27j. In fact, for a multiband system, 
the results of electrons tunneling and crossing the interface will be very direction dependent 
in terms of multiband FS topology. Moreover the existence of the surface bound states is 
related directly to the electron pairing potentials of different bands. Therefore a generalized 
theory to specifically study the possible midgap states in iron-pnictide superconductors and 
comparing to the PCARS measurements is in great demand. This is indeed the goal of 
the current paper. Here we will first show that zero-energy surface bound states exist in 
the possible s^-wave Fe-pnictide superconductors. While for other pairing symmetries such 
as (i-wave that support zero-energy bound states, the window of interface orientation for 
observing these states is usually large 16[, here in contrast, we find that zero-energy bound 
states can exist only when the orientations of interfaces fall into a small window around ±6° 
near the [100] direction. Furthermore, off the [100] direction, the zero-bias peak disappears 



and is replaced by two coherent peaks, which is due to directional dependence of QPs 
interplaying between different bands. These features provide unambiguous signals to test 
the possible s^-wave pairing in Fe-pnictide superconductors. 

The Model — For Fe-pnictide superconductors, the so-called ai and a2 Fermi sheets are 
concentric and nearly circular hole pockets around the F point. While /3i and ^2 Fermi 
sheets are nearly circular electron pockets around the M points 
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30| . High-resolution 



angle-resolved photoelectron spectroscopy (ARPES) measurements have been performed on 
superconducting Baa;Ki_2,.Fe2As2 from under-, optimally- to over-doped samples. It was 
shown that SC gaps on each FS are nearly isotropic. Besides, a large gap was observed on 
the two small hole-like /3i and FSs as well as on electron-like ai Fermi sheet, while a 
small gap was observed on the large hole-like ^2 FS. These FS sheets are sketched in Fig. [H 
Our discussions will follow the above picture of Fermi surfaces and SC gap structures. 

We first consider a Fe-pnictide bulk superconductor with a perfectly fiat and infinitely 
large interface located at x = 0. For each band, QP states have a coupled electron- hole 



character and can be described by the BdG equations (31 1 
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where ^ = — /i^V^/2m — /i with /i the chemical potential and m the electron mass. We shall 
seek a solution of the form 



u (r) = e^^^^'ri (r) and v (r) = e'^^'^'x (r) 



(2) 



where ki? is Fermi wavevector and in contrast to the plane- wave exponential factors, rj (r) 
and X (i") s-re slowly varying functions. Neglecting the second derivatives, we thus obtain 
the Andreev equations 

[ih {wf ■ V) + £;] 7/ (r) + a (r) x (r) = 0, 

[zn(v^- V)-E]x(r) + A(r)r7(r) = 0. (3) 

It is assumed that the pairing gap function takes the steplike form, A(x, = Ae^'^9{x). 
In reality, it sags a little near the interface at the distance of about the mean free path. In 
the following discussion, we neglect the above so-called proximity effect since it brings only 
small corrections. As the wave- vector components, ky, parallel to the interface are conserved 
for all possible processes, the problem is effectively reduced to one- dimensional one. 
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FIG. 1. (Color online) Panel (a): schematic plot of the Saint- James cycle for an Fe-pnictide two- 
band superconductor interface. The interface is oriented along the [100] direction. Panel (b): 
schematic plot of the s-wave SC gap amplitudes associated with different Fermi sheets. 



The two bands are assumed to be decoupled completely and thus one can proceed to 
obtain various reflection coefficients for each individual band. Outside of the interface, 
electrons and holes, if existing, are free electrons and free holes and we denote their wave 
functions by with ky suppressed. Here e (h) denotes for electron (hole) and + (— ) 

corresponds to movement parallel (antiparallel) to the x-axis. Inside of the superconductor, 
one can solve the Andreev equation ([3]) to obtain electron- and hole-like quasiparticle (QP) 
wave function \l/^^^^(a;). 

When a free electron is injected into the interface, two possible reflections and their related 
coefficients are as follows: (a) normal reflection (reflected as electrons) with the coefficient 
Tee and (b) Andreev reflection (reflected as holes, due to electron and hole coupling in the 
k subspace) with the coefficient Vgh. On the contrary, when a free hole is injected into 
the interface from the vacuum, there also exist two kinds of reflection: (a) normal reflection 
with the coefficient Vhh and (b) Andreev reflection (reflected as electrons) with the coefficient 
The- Matching the wave functions and their derivatives at the interface (x = 0), after some 
tedious but straightforward calculations, the Andreev reflection coefficients can be obtained 
to be 

g_iarccosf ^ for ^ < A, 
g-arccoshf ^ > A. 



feh{he) 



(4) 



For the subgap (or midgag) state, < A, we thus have for the total Andreev reflection 



eh(he){E)\'^ = 1. The normal reflection can be generally neglected 
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32|. 



With various reflection coeflicients, we are in the position to discuss the QP tunnehng. 



Following the formulation of Hu 12|, surface bound states of a Fe-pnictide superconductor 
can be obtained by considering QP states of a metal layer of width w attached to the 
interface in the limit of — )■ 0. In this limit, the solution to the Andreev equations Q 
inside the metal layer for a given ky can be written as 

V'(0,E)=aV;+(0,E) + 6^A+ (0,E) 

+cC(0,i?) + #,"(0,i?), (5) 

Due to non- conservation of momentum perpendicular to the interface at x = when QPs 
scatter off the interface, QP will scatter across a and /3 bands as shown in Fig. [1] for the 
case of /Cy = 0: An electron in the metal layer moving towards the interface will be Andreev 
reflected as a hole by the pair potential in /3 band. This hole will then be specularly reflected 
at the outer (free) surface, after which it will be Andreev reflected as an electron by the 
pair potential in a band. This electron will then, in turn, be specularly reflected on the 
free surface, which will form a close Saint- James cycle [l^. Quantitatively, one decomposes 
wave functions into partial waves in a and (3 bands by writing "^A^/j) = i^aeih) + '^%{h) with 
the following relations 

cCe(0,i?)=r,erfCJO,i?), 

6V'JjO,E) = re,aV;;,(0,E), 
(iC,(0,E) = -6^+ (0,E), 

aV';,(0,E) = -cCe(0,i?), (6) 

where the last two equations result from the specular reflections on the free surface. Note 
that the reversed cycle is an equivalent cycle in which one simply exchanges a and /3 in (E]). 
Clearly, Eq. (|6]) implies r-e/i(/3)r/ie(a) = 1 and reh{o:)rhe{f3) = 1 for the reversed cycle. By 
using Eq. (|4]), we flnd that the midgap energy in the two-band system must satisfy 

En En 

arccos + arccos — ^ = ± (</)q, — (pis) + 2mi, 

Aa Ap 

with n = 0,±l,--- (7) 

Eq. (I7j) represents one of the major results in this paper. If one identifles the scattering 
phase across the interface as the generalized momentum, pi = ±0j + arccos Eq. (HI can be 
rewritten as Yli={a p) Pi ~ 2?^7r, which is in consistence with the semiclassical quantization 
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FIG. 2. (Color online) Panel (a)-(c): schematic plots of different-direction electron tunneling on 
iron-pnictide SC interfaces. The arrows denote the electron tunneling directions and how different 
bands are coupled. The arrows with red sold lines show the particular directions where zero-energy 
Saint-James cycle can form. The weights of the arrows denote their possible contributions to the 
differential conductance. Panel (d)-(f): fitting to the G-V data of (d) SmFeAsOo.gFo.i sample 



(Tc = 51.5K) (taken from Ref. 



23|]); (e) TbFeAsOo.gFo.i sample 



and (f) SmFeAsOo.gFo.i sample (T^ = 51. 5K) (taken from Ref. [23]) 



m = 50K) (taken from Ref. 



24); 



condition. When Acf) = (j)a — (t^p = ±7r and En < Aq, < A/3, Eq. ((Tj) only supports the 
zero-energy solution {En = 0). It indicates that zero-energy surface bound states exist in 
the semiclassical approximation. Furthermore, since both Saint- James cycle and its reversed 
cycle support zero-energy state, there will be two non-dissipative currents associated with the 
zero-energy states. It should be noted that no Saint- James cycle exists when Aq, < En < Ap 
or Aq < A^ < En because according to Eq. (jlj), the currents will decay exponentially. 

Results and Discussions - To consider how the above midgap states are related to real 
PCARS measurements, we extend the above analysis to cases with /cy 7^ in the following. 
If the pairing gaps between different FSs are extended s-wave with sign reversal, there 
will be four kinds of zero-energy surface bound states (or zero-energy Saint- James cycles) 
corresponding to various combinations of interband FSs: ai — /3i, ai — /32, «2 — and 
02 — ^2- In addition to the above four zero-energy surface bound states, there also exist 
nonzero-energy surface bound states (or nonzero-energy Saint- James cycles) corresponding 
to various combinations of intraband FSs: ai — ai, 02 — f^i — f^i, and (32 — /32- 

To extend our analysis to investigate PCARS, the metal layer is replaced by a metal tip 
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occupying x < 0. When an electron is injected into the interface from the tip, there occur 
two possible reflections: (i) normal reflection (reflected as electrons) with the coefficient 
and (ii) Andreev reflection (reflected as holes) with the coefficient r^- In terms of r^r and r^, 
the resulting wave function for the normal side (x < 0) can be written as \t'Ar = ip^ +ri\]il)~ + 
f'A'iph- side, for the zero-energy Saint- James cycle, the transmitted electron-like 

and hole-like QPs are mainly in different bands, experiencing different pairing potentials 
for s^-wave Fe-pnictide superconductors. We shall take the electron-like QP to be in the a 
bands. The corresponding wave function is obtained by solving the Andreev equation 
([3]) with pairing potential A^e*'^". Similarly, the hole-like QP wave function \t'^ can be 
obtained by solving the Andreev equation ([3D with pairing potential A^e*'^''. Superposition 
of the two eigenstates, \E'^ and \E'^, will give a resulting wave function "^s = ci'^t + 
for the SC side. Matching the wave functions and their derivatives at the interface x = 0, 
he normal reflection amplitude rjv and the Andreev reflection amplitude ta can be found 
1^. Similar calculations apply to non-zero Saint- James cycles. From the knowledge of r^r 
and ta, the normalized differential conductance dl/dV is found according to the formula 

as = l + \rA? -\rN?m- 



When the interface is along the [100] direction as shown in Fig. [T] or Fig. [2](a), the 
dominant contribution to the current comes from electrons with ky ~ 0, i.e., most electrons 
tunnel normal to the interface. Due to non-conservation of kx-, QPs will scatter into different 
bands. The low-energy differential conductance will thus be dominated by the zero-energy 
bound states arising from the interband QP transitions although the nonzero-energy bound 
states may also have some effects on it. Thus the ZBCP found in PCARS is due to the 
zero-energy bound states that exist in the interface when electrons are tunneled normal to 
the interface along the [100] direction. This is verified in Fig. [2]^d), where we have simulated 



the PCARS data presented in Ref. 23| by assuming that all electrons are tunneling across 
the interface along the [100]-direction. The best fitting is obtained with Aq, = 8 and A^ = 8 
under the proposed sign-reversal s^-wave pairing. The normalized barrier height is taken to 
he Z = 5 and the scattering broadening P = 0.5. It should be emphasized that the ZBCP is 
more sensitive to the phase difference than the gap amplitude and when a more realistic FS 
together with a more realistic extended s-wave gap are used, a much better fitting, especially 
in the nonzero-energy region, will be obtained. The ZBCP found here can only exist when 
range of k„ covers 



window of ±6° for interfaces around [100] direction. 

We next consider the case when the interface is off the window of ±6° about [100]. In 
Fig. Et^b), we consider tunnehng in the [110] direction with corresponding range of ky. The 
directional detuning results in the absence of ZBCP. However, we found that the differen- 
tial conductance is dominated by the nonzero-energy bound states due to intraband QP 
transitions. This is the key why two coherent peaks are often observed in the PCARS mea- 
surements. In Fig. El^e), we use the two-gap s^'^-wave model to fit the PCARS data reported 



in Ref. [2J]. The best fitting was obtained with the two gap amplitudes set to = 8.5 and 
= 5. The broadening is taken to be F = 0.5. 

In Fig. [2](c), we consider the case when the orientation of the interface is between [100] 
and [110] but falls into the window of ±6° about [100]. In such case, both nonzero- and 
zero-energy bound states are equally important to the differential conductances and ZBCP 
is often observed due to zero-energy midgap state. Moreover, when electrons tunnel into the 
superconductor across the interface, QPs from different bands are coupled. Due to different 
gap amplitudes on different FSs and finite QP life time, one effective gap amplitude can 



be generally observed in experiment. In Fig. [2](f), we 
reversal [as done in Fig. [2](a)] to fit the PCARS data 



lave used one gap amplitude sign 
For the best fitting, the effective 



gap amplitude is taken to be A = 6.5 with broadening being F = 1.5, while the sign-change 
gap amplitude is taken to be A = 5 with the broadening being F = 0.01. 

In summary, we have investigated signatures of antiphase s^'^-wave pairing in Fe-pnictide 
superconductors due to zero-energy surface bound states. It is found that while for other 
pairing symmetries such as rf-wave that support zero-energy bound states, the window of 
interface orientation for observing these states is usually large, for s^'^-wave, however, zero- 
energy bound states can exist only when the orientations of interfaces fall into a small 
window around ±6° near the [100] direction. Off the [100] direction, the zero-bias peak 
disappears and is replaced by two coherent peaks, due to directional dependence of QPs 
interplaying between different bands. Our results give a unified explanation to the PCARS 
experimental data in various directions and indicate strongly that the pairing of Fe-pnictide 
superconductors is the s^'^-wave in different bands. 
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